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Abstract 



A method for generating irreducible polynomials of degree n over the finite field 
GF(2) is proposed. The irreducible polynomials are found by solving a system of 
equations that brings the information on the internal properties of the splitting field 
GF(2 n ) . Also, the choice of a primitive normal basis allows us to build up a natural 
representation of GF(2 n ) in the space of n-binary sequences. Illustrative examples 
are given for the lowest orders. 
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1 Introduction 



The problem of determining the irreducible polynomials over a finite field of order 
q, GF(q), is a task, from a qualitative mathematical point of view, similar to that 
of finding the prime numbers in the set of integer numbers [1, 2]. Such polynomials 
are necessary in numerous applications. For instance, the generation of algebraic 
error-correcting codes and cryptography [3] or the construction of fast transforms are 
based on finite field properties [4] that require algorithms to obtain the irreducible 
polynomials f n (x) of degree n over a particular finite field. 

Tables of irreducible polynomials are found in the literature [5]. These can be 
built up by different methods [6, 7], including a close analog to the 'sieve' method 
used in making tables of prime numbers. In this work a method to generate the 
irreducible polynomials of degree n over the field GF(2) is proposed. The strategy 
comprises the construction of a representation of the splitting field GF{2 n ) over the 
set of binary n-digits sequences. Each periodic orbit and its m permutations, with 
m a divisor of n, identify a irreducible polynomial of degree m and its m roots on 
GF(2 n ). 

In section // the connection between the field GF(2 n ) and the space of n-binary 
sequences, B n , is identified. The two internal operations over B n necessary in order 
to obtain a particular representation of GF(2 n ) are defined in section 111. The 
method to find the irreducible polynomials f n (x) over GF(2) is disclosed in section 
IV. Section V summarizes the general scheme to be followed in order to obtain 
f n {x) for all n. Two illustrative examples are given in section VL Last section is 
concerned with our conclusions. 



2 GF(2 n ) and B n 

Notation: Let us call B n the space of binary sequences of n digits. B n has 2 n 
elements. Each element O G B n is said to be an orbit. If O is given by the binary 
sequence (a n _io; n _2 • • • (Xiao), = or 1), the left shifted orbit O' defined 
by the sequence (a ra _2«rt-3 * ' ' ai«o a ?i-i) is represented by shift(0). If m is the 
smallest integer such that O = shift m (0) = shift ( — ( shift ( O) ) ) then the 
orbit O is said to be of period m or m-periodic, and it is represented by O m . It is 
obvious that if O m G B n then m is a divisor of n, m | n. The set [O m ] formed by 
the m different shifts of an m-periodic orbit, O m , is called an orbital of period m: 
{O m } = {O m ,shift(O m ),shift 2 (O m ) r ■ ■ ,shift m ~ l {O m )}. As an example, [0 3 ] = 
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[100] = {100,010,001} or [0 2 ] = [0101] = {0101,1010} = [01] = {01,10}. 



Fundamental Idea: Firstly, let us note that the number N n of irreducible 
polynomials (LP.) of degree n over GF(2), is equal to the number of orbitals [O n ] 
of period n (table 1,[8, 9]): 

= 2 " ~ £ ' =1 m ' Nm , (1) 

n 

where {mi, rri2, . . . , rrik} are the integer divisors of n less than n. 
The idea is to build up a representation of GF(2 n ) over B n , with the constraint that 
each orbital [O m ] contained in B n has an associated LP. , f m (x), of degree m over 
GF(2). Moreover, the m orbits O m that expand the orbital [O m ] are the m roots of 
fm( x ) hi GF(2 n ). Thus, two internal operations (+, *) are defined over B n in such a 
way that (£>„, +, *) be the splitting field for all LP. f m (x) over GF(2) , with m \ n. 



3 Field Construction in B 



Choice of a Basis: If we consider GF(2 n ) a vectorial space over GF(2), a primitive 
normal basis, B an ={a Q , a±,..., o" n _i}, can be found in GF{2 n ) [10, 11]. B an to be 
normal means that of = mod «• B an to be primitive means that ao generate 

GF(2 n )* = GF{2 n ) — {0}. Therefore, if we denote by (+, *) the two internal opera- 
tions in GF(2 n ), each element a G GF(2 n )* has a unique representation as addition 
(+) or as product (*) of the B an elements (tables 2, 3, 4, 5, 6): 

(i) + -representation: each element a G GF(2 n ) can be expanded in a linear com- 
bination of B na , a = J27=o a i a i with ctj G GF(2) = {0,1}. The binary sequence 
a + = (a n _io; n _2 • • -cnao) £ £>n is the +-representation of a. 

(ii) * -representation: each element a G GF(2 n )* is generated by a , consequently 
a = <7q with G {1,2, •••,2™ — 1}. If we take the binary expression of the ex- 
ponent, k = YJlZo (3i2\ the binary sequence a* = (/5 n _i/3 n _ 2 • • • PiPo) G B n is the 
*-representation of a. It can also be reformulated as a = n^r Vf where a? = 1. It 
is obvious that (a 2 )* = shift(a*). 

Example (table 5): If -B -4={o r o ) 0i, a"2, (T3} is a primitive normal basis, the 
following expansions are obtained for some elements of GF{2 A ): 

a = <7q = o"o + 0i + (73 — > a + = 1011, 
a = ctq = o"o * 0i - ► a* = 0011, 



6 = a % = 02 -»• &+ = 0010, 
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b = a\ = a 2 -> K = 0010, 



c = a\ = a 2 + 0-3 — > c + = 1100, 
C = (Tq = (To * en * °"2 — ► = 0111. 

The particular elements, a = = a% and a = Gq _1 = cr n _i * o" n _ 2 * • • • * o\ * Oq — 
&n-i + °"n-2 + • • • + o"x + o" = 1, verify a + = a*. In this case, (00 • • • 010 • • • 00) and 
(11 • ■ • 11 • • • 11), respectively. 

Internal operations in B n : we explain now how to perform the two internal 
operations (+, *) in B n . 

+- operation: the property o"j + o"j = (the characteristic of GF(2 n ) is 2) imposes 
the algorithm for performing this operation. Given a, b G GF(2 n ) the element a + b 
is obtained by summing bit to bit the binary sequences of a and b (mod2) in the 
+-representation. If we have 

a + = a n -ia n -2 ■ ■ ■ a\a , 

b+ = Pn-lPn-2 ' ' ' PlPo, 

(a + b) + = 7„_i7„_2---7i7o, 
the addition + presents the form: 

(a + b) + = a + + b + mod 2, 
7i = (a* + A) mod 2. 

(/?„, +) is a representation of the abelian additive group defined over GF(2 n ), where 
the neutral element is = (00 • • • 00) and the inverse element a of a is the element 
itself, a + a = 0, then a — a. 
Example (table 5): 

(a 3 Q ) + 1011, 

(^ 4 )+ - 0010, 
(a 3 + a 4 ) + -> 1001 = <x 14 . 



*- operation: the property of = <7(i+i) mo dn gives the clue to perform this opera- 
tion. Given a, 6 e GF(2 n ) the element a*6 is obtained after performing the addition 
with carrier (mod2) of the binary sequences of a and b in the ^representation. If 
we have 

a* = a n -ia n -2 ■ ■ ■ a\a , 

b* = (5 n -\fl n -2 • • • P1P0, 

(a + 6)* = 7 n _i7„_ 2 • • • 7i7o, 
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then the multiplication * presents the form: 



(a * b)* = a* * = (a* + 5* with carrier) mod 2, 
7i = (aj + /3j + carrier^) mod 2. 

To start we suppose carrier - = 0. If a n _i+ / 9 n _i + carrier n _ 1 > 1 then the operation 
is restarted with carrier^ = 1. (£> n , *) is a representation of the abelian multiplica- 
tive group defined over GF(2 n )*, where the neutral element is 1 = (11 • ■ ■ 11) and 
the inverse element a -1 of a is represented by the complementary binary sequence, 
(a -1 )* = cp(a*) given by digiti^a' 1 )*) = (digiti(a*) + 1 mod 2). This is a conse- 
quence of the fact that if a = <Tq then a -1 = oq 2 ^ k . The multiplication by zero 
is (a * 0)* = a* * (00 • • • 00) = (00 • • • 00) = by definition. It can also be seen that 
if (a * ft)* = 0, then only one of the two sequences are null, a = or b = 0. 
Example: 

carrier^ = 0, 

((Tq 11 ), -> 1011, 

(a 4 )* -> 00J0, 

(<7 U *(7$), -> 1101 = (7 13 . 



carrier q = 1, 
(o-q 1 )* 
(1)* 

K 1 * i)* 



1011, 
1111 , 

1011 = (Jq 1 . 



a* = 1011, 
(a- 1 ^ = 0100. 

Connection between + anc? * : the central question in order to build up a represen- 
tation of the field GF(2 n ) over B n is to find the relationship between a + and a*. 
This construction can be performed if that relation is known for the elements (Tj * o~j 
because all the other elements can be reduced to products of <7j * Uj type. Thus, if 
we write 

n-1 

o-i*o-j = ^eijiai, 

1=0 

there are n 3 constants e^i (= 0, 1) defining the field GF(2 n ) completely. By con- 
struction <jj * G{ — o~i+i, then em = By the commutative property, = e^. 
Also, (<7j * <Tj) 2 = <7j + i * (7j + i, then shift((o~i * cr,-)*) = (cj+i * Cj+i)*- If we impose 
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this property for the +-representation, that is, shift((<7i * = (cj+i * + i 

then e( i+1 )(., +1 )(7 +1 ) = e^. Thus, only (n — l)/2 products <Ji * aj when n is odd, and 
n/2 products Oi * <jj when n is even, must be known in order to determine the field 
completely. That is, a total number of n 2 /2 or (n 2 — n)/2 constants e^, depending 
on the parity of n, defines the field. These constants cannot be freely chosen as the 
distributive property must be verified: 

a* {b + c) = a * b + a * c, 
(a* * (b+ + c + )*) + = (a* * 5*)+ + (a* * c*) + . 

All the constants stay fixed by these conditions in some cases and the natural 
construction (£>„, +, *) is unique. In general, this is not the case and there exist 
different constructions of GF(2 n ) in B n . 



4 Construction of the Irreducible Polynomials 

Let us recall that GF(2 n ) is the splitting field of all irreducible polynomials f m (x) 
over GF(2) of degree m < n and m \ n. If a e GF(2 n ) is a root of f m (x), the m roots 
of this irreducible polynomial are {a, 2m }. Therefore, the m orbits 

that integrate an orbital [O m ] G B n are the m roots of the irreducible polynomial 
f m (x) in the representation (B n ,+,*) of GF(2 n ), as explained in the last section. 
For instance, in the case n = 4 (table 5), the elements {ao&i, &1&2, C2C3, C3C0} are 
the roots of the irreducible polynomial x 4 +x 3 +x 2 +x+l. The objetive of this section 
is to establish a method to find the irreducible polynomials over GF(2) working out 
this property. 

Trace of an orbital: The trace of an orbital [O m ], Tr(O m ), is defined as the 
trace of one of its elements: {a, a 2 , a 22 , • • • , a 2(m x) } with a e GF(2 n ), that is, 

Tr(Oj=a + a 2 + a 22 + ... + a 2(m_1) . 

The trace is an element of GF(2), in this case or 1. The number of traces A^ ra ' ) in 
GF{2 n ) is equal to the number of orbitals [O m ] in i3 n , 

iV^ = ^2 N m with m \ n. 

m 

Let us call Nj?i the number of traces whose value is 1 and N^o the number of those 
whose value is 0. Hence, N^ n) = N$ + N$ . 

The important point for our proposal is that the coefficients of an LP. over GF(2) can 
be obtained by adding (mod2) the traces of different orbitals of the field GF(2 n ). 
Therefore the LP. can be identified if the traces are previously calculated. 
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Example (table 5): (i) the LP. associated to the orbital {<7o<7i, o~\o~ 2 , o~ 2 o~3, 0300} 
is a polynomial of degree 4, x A + cfe 3 + ex 2 + fx + 1. Its coefficients are functions of 
the traces: 

d= <Jq(Ji + o x o 2 + 0-20-3 + a 3 (T = TV(<7 Ci), 

e= f7 (Ji + (Ji(T2 + (7 2 (73 + (7 3 (To + (7 (Ji (72(73 + (7o(7lf72f7 3 = TV(<7 (7i) , 
/ = CToCl + 0"lO"2 + 0" 2 (7 3 + (7 3 (7o = TV(<7 <7i). 

(ii) the LP. associated to the orbital {a , <7i, a 2 , <r 3 } is also of degree 4, x 4 + ax 3 + 
bx 2 + cx + 1, where: 

a= o- + o x + cr 2 + 0-3 =Tr(a ), 

b= (7 (7i + 0-10-2 + cr 2 0-3 + 0-30-0 + a a 2 + (7i 0-3 = Tr(ao<7i) + Tr(a a 2 ) : 

f= (7 (7i(72 + (7i(7 2 0"3 + O"2O"3O"0 + O"3O"0O"l = TV(<7o<7i<7 2 ) . 

Newton's Formulas: the fundamental theorem on symmetrical polynomials can 
guide us in finding the different traces in the field GF(2 n ). This theorem can be 
expressed as follows [12]: Let r±, r 2 , ■ ■ •, r n be different variables and let us call 8\, 
■ ■ •, 8 n the following symmetrical polynomials in those variables: 

1 = n + r 2 H h r n , 

02 = nr 2 + rir 3 H h r n _ir„, 

3 = rir 2 r 3 + rir 2 r 4 H h r„„2^n-i^n, 

6» n = nr 2 r 3 ---r n . 

Then every symmetrical polynomial in the variables {ri, r 2 , • • •, r n } and with co- 
efficients over a field F can be expressed as a polynomial in the new variables 
{9i, 9 2 , ■ ■ ■ , 6>„} with coefficients in F. 

In the particular case of the following symmetrical polynomials s^: 

s k = r\ + r\ + • • • + r k n k = 1, 2, 3, • • • 

the expansion in the variables is given by the Newton's formulas [13]: 

si-0i = 0, 
s 2 -sA + 20 2 = 0, 

s 3 - s 2 e ± + Sl e 2 + 30 3 = 0, (2) 

s k -s k - 1 1 + s k _ 2 2 ---- + (-l) k - 1 s 1 k _ 1 + (-l) k (k)0 k = 0, 
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where Oj = for j > n and the symbol (k) is interpreted as (k) — 1 + 1+ • • • +1. 

In our problem, for instance, if we take = <7j_i the expressions of s& and 8^ 
are functions of the GF(2 n ) traces. Then the Newton's formulas give us a system 
of equations for the different traces of GF(2 n ). By solving this system it is possible 
build up the LP. of degree n over GF(2). 

Finding the Traces: If the elements {a , a 1: • • • , <r n -i} £ GF(2 n ) are taken as 
the variables {ri, r 2 , • • •, r n } in the equations (2), a system of equations for the 
different traces is obtained. The total number of unknowns (traces) is equal to the 
total number of orbitals, \ in GF(2 n ) (eq. 4). By construction: 

Tr(0) = 0, 
Tr(a ) = 1, 
Tr(a ai ■ ■ -cr n _i) = 1. 
For large n the number of orbitals is of the same order as 

n 

Noting that rf" = of™ = cr, = r { we will have 2 n equations. Therefore, a system 
of 2 n equations with 2 n /n unknowns (many equations are interdependent) should 
be solved. Only those solutions verifying that the number of traces with value 1 
is equal to N^i are accepted. If one of these solutions is chosen, then the LP. are 
directly calculated. 

Note: Let us recall a property that allows us to generate the LP. by pairs. 
Two polynomials, f(x) and g(x), are reciprocal if they can be expressed in the form: 

f(x) = a n x n + a n -\x n ~ x H V a±x + a , 

g(x) = a x n + a x x n ~ x H V a n ^x + a n , 

or equivalently, g(x) = x n f(x~ 1 ). This means that the roots of g(x) are the inverse 
roots of f(x). Consequently, if f m (x) is an LP. of degree m, the reciprocal polynomial 
g m (x) is also an LP. of degree m. This property can be used as another useful tool 
to facilitate the finding of all irreducible polynomials. It is necessary to remark at 
this point that both polynomials can be the same, that is, f m (x) = g m (x). 

5 General Methodology 

Summarizing, the steps to follow for building up a natural representation of the field 
GF(2 n ) in B n and for finding all LP. of degree n over GF(2) are the following: 
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(1) All the independent n-binary sequences (orbitals [O m ] with m \ n) in B n 
are found (see appendix). Its total number, N^f \ is the total number of LP. over 
GF(2). The splitting field is GF(2 n ). Let us observe that each orbital [O m ] has 
associated an LP. f m (x) of degree m. The roots of f m (x) in (£>„,+,*) are the m 
roots given by the m orbits of the orbital [O m ] . It also is possible to evaluate N^i 
by calculating directly the trace for each orbit O m , 

Tr(O m ) = O m + shift(O m ) + shift 2 (O m ) + ■■■ + shift m -\O m ), 

where shift(O m ) is the left shift of the orbit O m . If #i(O m ) represents the number 
of 1 in the orbit O m multiplied by m/n then it can be deduced that 

Tr(O m ) = 1 if #i(O m ) is odd, 
Tr(O m ) = if #i(O m ) is even. 

For instance, 



#i(0011) 


= 2 • 


4 
4 


= 2 


Tr(OOll) 


= 0, 


#i(0111) 


= 3 • 


4 
4 


= 3 


=>• Tr(Olll) 


= 1, 


#i(0101) 


= 2 • 


2 
4 


= 1 


=^ Tr(0101) 


= 1, 


#l(llll) 


= 4- 


1 

4 


= 1 


=>• Tr(llll) 


= I. 



(2) To calculate the trace associated to each element in GF(2 n ), the 2 n non- 
independent equations deriving from the Newton's formulas are established. There 
are ~ 2 n /n unknowns. 

(3) A solution of the former equation system is worked out. This solution must be 
compatible with the constraints: (i) Tr(0) = 0, Tr(o"o) = 1 and Tr(aoa\ ■ ■ -o" n _i) = 
1; (ii) the number of traces with value 1 must be equal to (calculated in 
point (1)). The different possible solutions are related with the different possible 
constructions of GF{2 n ) in B n . 

(4) The irreducible polynomials are explicitely found in a deterministic way from 
a particular solution obtained at point (3). 

(5) The construction of (B n ,+,*) is performed from the LP. calculated in (4). 
This is a natural representation of GF(2 n ). 

(6) Note : As it can be seen in the tables below, an only representation of GF(2 n ) 
in B n is possible when {n = 1, 2, 3, 4}. 
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6 Examples 



In this section the process formerly explained is applied for the cases n — 4, 5. 

n—4: Let us suppose GF(2 4 ) expanded by the primitive normal basis B a4 = 
{(T , 0i, o- 2 , cr 3 }. Then of = <7( i+ i) mod 2 . The field's elements {0, cr , 0"q, <7q, • • • , a 1 , 4 , cr 1 , 5 = 
1} are generated by <jq. If an element a G GF(2 4 ) verifies a 2 ™ = a, the m pow- 
ers of a, {a, a 2 , a 2 , • • • , a 2 " 1 }, are the m roots of an LP. of degree m over GF{2). 
Obviously, m is a divisor of 4. The different representative elements of this kind in 
GF(2 4 ) are: 



a G GF(2 4 ) 


m 





1 




4 


a a 1 


4 




2 




4 




1 



The different orbitals [O m ] of period m, with m | 4, in i3 4 are given by the represen- 
tative orbits: 



[O m ] 


Tr(O m ) 


m 


0000 





1 


0001 


1 


4 


0011 





4 


0101 


1 


2 


0111 


1 


4 


1111 


1 


1 



Thus, Np = Ni + N 2 + N 3 = 6 and iVf/ = 4 (see Table 1). It can also be 
observed that the orbit {o- o" 2 , o^o^} is of period 2, then its binary +-representation 
is {0101, 1010} and by direct calculation Tr(<jQ<j 2 ) = 1- 

If {a , <7i, cr 2 , 0-3} are taken as the roots of an LP. of degree 4, the equations derived 
from Newton's formulas (where 9k = for k > 4) for the different traces can be 
established. The variables and are expressed as follows: 

9\ = (T + <7i + 02 + 0-3 

= 7V(<7 ) = 1, 
^2 = o-qO - ! + aiO- 2 + o- 2 o- 3 + a 3 o- + o- a 2 + o-xo- 3 
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= 7V(<7 ci) + Tr(a a 2 ), 

#3 = (7 (7i (T2 + (7i(72(73 + Cr 2 (73(7o + (73(7o(Ji 

= Tr(a aia 2 ), 
84 = a aia 2 as 

= Tr((7 (7i(72(73) = 1, 

Sk = cr + a\ + a\ + cr* 

= ^Tr(V) mii/i fc = 1,2, ■ --,15. 

i 

The independent Newton's formulas from the system (2) are: 

Tr(a a 2 ) = 1, 
Tr(a a 2 ) + Tr(a a 1 a 2 ) = 1, 
Tr((7 (7i) = 1. 

The solution is unique: 

Tr(a ) = 1, 

7V((7o<7i) = 1, 

Tr((x (T 2 ) = 1, 

7V(<7 <7i<7 2 ) = 0, 
Tr((7 (7i(72(73) = 1. 



Finally, the LP. are: 

[0 4 ] = {o" } -> + ax 3 + bx 2 + cx + 1 ^x 4 + x 3 + l, 

[O4] = {c^ 1 = 0Wi°2} — ► x 4 + cx 3 + 5x 2 + ax + 1 — > x 4 + x + 1, 

a = Tr(a ) = 1, 
6 = Tr((7 <7i) + Tr(o- o-2) = 0, 
c = Tr(cr (7i(7 2 ) = 0, 

[0 4 ] = {o- o"i} -»• x 4 + cfe 3 + ex 2 + ,/x + 1 ^x 4 + x 3 + x 2 + x + l, 

d = Tr((7 (Ji) = 1, 
e = Tr((7 (7i) = 1, 
/ = Tr((7 (7i) = 1, 

[0 2 ] = {<7 (7 2 } -»• x 2 + #x+l ^x 2 + x + l, 

5- = Tr(a o- 2 ) = 1, 
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[Oi] = {o- o"iO"2O-3} -> a; + 1 -> x + 1, 

[Ox] = {0} -> x -> x, 

where {ov--<7fc} is representing all the GF(2 4 ) elements that correspond to the 
orbital [O m ]. 

The step to find the +-representation of every a e GF(2 4 ) is the following. 
Let us take cr with associated LP.: x 4 + x 3 + 1, then: 

a* + a 3 + l = 0, 

a 4 = a 2 -> (<t{}) + = 1 + ((7 2 ) + = 1011. 

Taking the shifts of 1011, we find the +-representation for all the elements of the 
orbital {o"q} = {a ai}. The distributive property determines the +-representation 
of the orbit {oocrio^}, with the final result shown in table 5. 
Let us observe that this construction of GF(2 A ) is unique. 

n=5: We repeat the same process as one of the case n = 4. Let us suppose 
GF{2 b ) generated by the primitive normal basis B a5 = {a , a±, a 2 , cr 3 , 04}. Hence 
of = <7(j+i) m od 2- The field's elements {0, (Jo, <Tq, <Jq, ■ ■ ■ , erg , ^o 1 = 1} are generated 
by (Tq. If an element a G GF(2 5 ) verifies a 2 " 1 = a, then m \ 5 and the m powers of 
a, {a, a 2 , a 22 , • • • , a 2 ™ are the m roots of an LP. of degree m over GF{2). The 
representative different subsets of this kind in GF{2 b ) are: 



a E GF(2 5 ) 


m 





1 




5 


<7 <7i 


5 


O"0°"2 


5 


O"0°"l°"2 


5 


O"0°"l°"3 


5 




5 


O"0O"lO"2O"3O"4 


1 



The different orbitals [O m ] of period m, with m \ 5, in B$ are given by the represen- 
tative orbits: 
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\o 1 


Tr(0 \ 


1 VI 

ill 


(J (J (JUL) 





1 
1 


UUUU1 


1 





nnm 1 


n 

u 


O 


00101 





5 


00111 


1 


5 


01011 


1 


5 


01111 





5 


11111 


1 


1 



Thus, N£' = Ni + N 5 = 8 and Nfi = 4. 

If {<7 , &i, &2, 03, 04} are taken as the roots of an LP. of degree 5 the equations 
derived from Newton's formulas (where 9k = for k > 5) for the different traces can 
be established. The variables 9 k and Sk are expressed as follows: 

9\ = o-q + ai + a 2 + <r 3 
= 7V(<7 ) = 1, 

92 = 0"oO"l + 0"1<T2 + CT2C"3 + °"3°"4 + O"4C"0 

+Cr a 2 + (Ti(T3 + 0" 2 (T4 + <7 3 <7 + (74<7i 

= Tr(cr cri) + Tr(a <7 2 ), 

93 = (7o(7lf72 + (Tl(J2(J3 + (72(73(74 + C3(T4(Jo + (J4(To(Ti 

+O- O- 1 O- 3 + (7i(72(74 + 0- 2 <7 3 (7o + <7 3 <74<7i + er 4 er er2 

= Tr (a Q ai a 2 ) + Tr (owi^), 

#4 = O'0°'l o "2O"3 + (7l<7 2 (73(74 + 0"2<7 3 <740"o 
+O" 3 O"4O" O" 1 + 0"4<7o<7i<72 

= Tr((7o(7i 0-2(73), 

6» 5 = Tr(<7o<7i<7 2 <7 3 (74) = 1, 

Sfc = a fc + of + (7 2 fc + o| + a* 

= ^Tr(V) luit/i A; = 1,2, •••,31. 

i 

The independent Newton's formulas from the system (2) are: 

Tr(a a 2 ) +Tr(a a 1 a 2 ) +Tr(a a 1 a 3 ) = 1, 
7V(<7o<7i<7 3 ) + Tr(o- o-io- 2 o-3) = 0, 

7V(<7 <7i) + 7V(<7 <7i<7 2 ) = 0, 

Tr(a- a- 2 )(l+Tr (<7 <7i<7 3 )) = 0, 
Tr(<To<T 2 )(l + Tr((To<Ti)) = 0. 
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One solution compatible with the constraints is: 



1 r l°oJ 


— 


2V(<7 <7i) 


= o, 


Tr((7o(7 2 ) 


= o, 


Tr((T f7i(72) 


= o, 


Tr(cr cri(73) 


= 1, 


Tr(o- o-io- 2 o-3) 


= 1, 


Tr(o- o-io- 2 o- 3 o-4) 


= I. 



Finally, the LP. are calculated from this solution: 

[0 5 ] = {a } -> x 5 + ax 4 + 6a; 3 + cx 2 + dx + 1 -> x 5 + x 4 + x 2 

+ 1, 

[0 5 ] = {o^ 1 = a-oo^o^o^} — > x 5 + <ix 4 + cx 3 + bx 2 + ax + 1 — > x 5 + x 4 + x 3 

+X+ 1, 

a = Tr(o- ) = 1, 
b = Tr(a a 1 ) + Tr((io(72) = 0, 
c = Tr(a aia 2 ) + Tr(<To<7i<7 3 ) = 1, 
d = Tr(o- o-io- 2 o- 3 ) = 1, 

[0 5 ] = {o- o-i} -> x 5 + ex 4 + fx 3 + gx 2 + hx + 1 ^x 5 + x 3 + l, 
[0 5 ] = {(oWi) -1 = 0-00-10-2} -> x 5 + /ix 4 + #x 3 + /x 2 + ex + 1 ^x 5 + x 2 + l, 

e = 7V(<7 <7i) = 0, 
/ = Tr(croO- 2 ) + Tr(<7 O Ci0"20"3) = 1, 
5 = Tr(o- o-io- 3 ) + Tr(o- ) = 1, 
ft, = Tr(a a 1 a 2 ) = 0, 



[O5] = {0W2} 

[O5] = {(o- o- 2 ) _1 = 0-00-10-3} 

3 = Tr(a a 2 ) = 0, 
k = Tr(a ai(T2cr3) + Tr(o- o-io- 2 ) = 1, 
I = Tr(a ) + Tr(a ai) = 1, 



x 5 + jx 4 + /cx 3 + lx 2 + mx + 1 — > x 5 + x 3 + x 2 

+x + 1, 

x 5 + mx 4 + /x 3 + A;x 2 + jx + 1 -> x 5 + x 4 + x 3 



, 2 , 1 
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m = Tr(a cr 1 i7 3 ,) = 1, 

[Oi] = {(To (71(72 (T3 (74} — > X + 1 — > X + 1, 

[Oi] = {0} -> X -> X, 

where Wi---cFk} is representing all the GF(2 5 ) elements that correspond to the 
orbital [O m ]. 

The steps to find the +-representation of every a G GF(2 5 ) are the following. 
Let us take <7 with associated LP.: x 5 + x 4 + x 2 + x + 1, then: 

(Jo 5 + o£ + a\ + a + 1 = 0, 

a 2 = (n ; <r$ = (7 2 -> ((750+ = 1 + (cr ) + + (a 1 ) + + (<7 2 )+ = 11000. 

The shifts of 11000 are the +-representation for all the elements of the orbital 
Wo) = Wo^}- 

If we now take Uq with associated LP.: x 5 + x 3 + 1, then: 

K 3 ) 5 + K 3 ) 3 + i = 0, 

((7 9 ) + = 00110 - ((7 15 ) + = 1 + ((7 9 ) + = 11001. 

The shifts of 11001 are the +-representation for all the elements of the orbital 
{°o 5 } = Wocricr2cr 3 }. 

Proceeding now with Oq that has an associated LP.: x 5 + x 2 + 1, then: 

K 7 ) 5 + K 7 ) 2 + i = 0, 

°f = 4 - K 14 )+ = i + K 4 )+ = 11011. 

The shifts of 11011 are the +-representation for all the elements of the orbital 
{°o} = Wocricr 2 }. 

Taking <Tq 5 with associated LP.: x 5 + x 4 + x 3 + x + 1, then: 

(0 5 + (0 4 + (0 3 + ^ 5 + i = 0, 

4° = 00111 ; = 11011 - ((7 75 ) + = ((7 13 ) + = 

i + (^ 5 ) + + K 45 ) + + K 60 ) + = 11010. 

The shifts of 11010 are the +-representation for all the elements of the orbital 
Wo 1 } = {o"oO"iO"3}- 

Finally, the +-representation for the elements of the orbital = {o"o(7i} are the 
shifts of 01001. Other solutions for the Newton's formulas can allow different natural 
representations of the field GF(2 5 ) in B n . Table 6 is one of them. 
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7 Conclusion 



A possible method to build up the irreducible polynomials of degree n over the field 
GF(2) has been exposed in this paper. The strategy exploits the properties and the 
connection between the splitting field GF(2 n ) and the space of n-binary sequences 

The choice of a primitive normal basis allows us to perform a natural construction 
of GF(2 n ) in B n : (i) Two natural operations (+, *) are defined on B n in order to 
obtain a representation of the splitting field GF(2 n ); (ii) A system of equations 
for the different traces of the field is derived from the fundamental theorem on 
symmetrical polynomials; (iii) Different constraints, that have been established, limit 
the number of possible solutions of these equations; (iv) Any of these solutions gives 
us the irreducible polynomials of degree m, with m \ n; (v) The field (£>„, +, *) is 
finally built up with the help of the irreducible polynomials. 

We want to remark the direct relationship that seems to exist between the irre- 
ducible polynomials over GF{2) and the periodic binary orbits, in such a way that 
this is a one-to-one application. Following this ideas and under a similar scheme, 
this method could be applied to find the irreducible polynomials of degree n over 
a general finite field GF(q) by using the connection between the properties of the 
splitting field GF(q n ) and the space of n-q-ary sequences. 

Acknowledgements: This work was performed during a postdoctoral stage at 
the Ecole Normale Superieure of Paris in 1995. I thank J-M. Couveignes (Paris) and 
J. Cabeza (Zaragoza) for useful discussions and references on this subject of finite 
fields, and the Gouvernement Frangais for a research grant. 

Note : This manuscript should be understood as the approach of a 'newbie', plenty 
of questions, to a subject very far from its common working topics. 
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Appendix: Calculation of Nj, and N : 



(n) 



Tl 



An algorithm to calculate Nj? and iVf| J is proposed in this appendix. The 
process to find the different orbitals in £>„ involves the following steps: 

(1) The decimal value of the different n-digit orbitals is stored in an indexed 
integer variable 0(i), then < i < Nj?\ The initial conditions are: O(0) = 0, 
0(1) = 1. Another integer variable O p is defined and initialized as O p = 1. 

(2) The operational loop to find 0(k + 1) when all orbitals 0(m) are known for 
m < k is the following: 



(**) O p = O p + 1 
do i = 0, n — 1 

Op = O p ■ T mod (2" - 1) 
do j — 1, k 

if (Op~.eq.O(j)) goto (**) 
end do 
end do 

0(k + l) = O p 
k = k + 1 
goto (**) 



(3) This is repeated for all sequences of B n while O p < (2 n — 1). The maximum 
index k from the step (2) is k max = — 1. 

(4) The period (m) of each orbital 0(i) is calculated: 



per(0(i)) = m <^ m is the smallest integer 1 < m < n — 1 



(5) Each orbital is converted to its binary expression and its trace is calculated: 



Tr(0(i)) = 7V(0 6 (i)) = O ft (i) + shift(O b (i)) + shift 2 (O b (i)) + ■■■ + shift m (O b (i)) 



(6) The number of traces equal to 1 when i ranges from 1 to Nf> is N : 



0(i) = 0(i) ■ 2 m mod (2 n - 1) 




binary (O(i)) 
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Table Captions 



Table 1. Functions over the natural numbers J\f, that appear in this work, are 
calculated for the first n e N. N n is the number of irreducible polynomials of degree 
n over GF(2). Nj^ = J2 m N m with m\n represents the total number of orbitals in 
B n . N$ and N$ are the number of orbitals in B n with trace 1 or 0, respectively. 

Table 2,3,4,5,6. Natural representation of the finite field GF(2 n ) for the cases 
n = 1,2,3,4,5 in the space of binary sequences B n , where the two internal oper- 
ations are performed in the following way: (i) addition: (a + + b + ) + = a + + b + ; 
(ii) multiplication: (a * 6)* = a* * 6*. The set (B n ,+,*) is a representation of the 
finite field GF(2 n ). The m orbits that compose each orbital [O m ], with m \ n, are 
the m roots of an irreducible polynomial (LP.) of degree m over GF{2). 
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n 






JV T1 


/V (riJ 

iv T0 


1 


2 


2 


1 


1 


2 


1 


3 


2 


1 


3 


2 


4 


2 


2 


4 


3 


6 


4 


2 


5 


6 


8 


4 


4 


6 


9 


14 


8 


6 


7 


18 


20 


10 


10 



Table 1: 



GF(2) 


element: a 


a* 


a + 


Trace 


LP. 


c 




1 


1 


1 


x + 1 




















Table 2: 



GF(2 2 ) 




element: 


a 


a* 


a + 


Trace 


LP. 


c 


c 


c 


01 


01 


1 


x 2 + x + 1 


c 2 


Cl 


Cl 


10 


10 






c 3 


c + Cl 


<7 (7i 


11 


11 


1 


x + 1 











00 


00 








Table 3: 
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GF(2 3 ) 


element: a 


a* 


a + 


Trace 


LP. 


00 
< 


00 
01 


(Tl 
02 
03 


001 
010 
100 


001 
010 
100 


1 


x 3 + x 2 + 1 


-o 6 
O 5 


O"0 + 0" 2 
O"! + (T 
a 2 + (Ti 


0001 
0102 
(T2(To 


Oil 
110 
101 


101 
Oil 
110 





X 3 + X + 1 


-o 7 


(To + (Tl + 0"2 


000102 


111 


111 


1 


X+l 











000 


000 





X 



Table 4: 



GF(2 4 ) 


element: a 


a* 


a + 


Trace 


LP. 


(To 


00 


00 


0001 


0001 






0o 2 


01 


01 


0010 


0010 


1 


x 4 + x 3 + 1 


0o 4 


02 


02 


0100 


0100 






*o 8 


03 


03 


1000 


1000 






0o 3 


00 + 01 + 03 


0001 


0011 


1011 






-o 6 


01+02 + 00 


0102 


0110 


0111 


1 


x 4 + X 3 + x 2 + X + 1 




02 + 03 + 01 


0203 


1100 


1110 






0o 9 


03 + 00 + 02 


0300 


1001 


1101 






0o 5 


00 + 02 


0002 


0101 


0101 


1 


x 2 + X + 1 




01 + 03 


0103 


1010 


1010 






0o 7 


02 + 03 


000102 


0111 


1100 






< 


03 + 00 


010203 


1110 


1001 





x 4 + X + 1 




00 + 01 


020300 


1101 


0011 






0O 11 


01 + 02 


030001 


1011 


0110 






CT 


00 + 01 + 02 + 03 


00010203 


1111 


1111 


1 


X+l 











0000 


0000 





X 



Table 5: 
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GF(2 5 ) 


element: a 


a* 


a 4. 


Trace 


LP. 


(Tn 
u U 


(Tn 
u U 


(Tn 
u (J 


00001 


00001 











<Ti 


(Tn 


00010 


00010 




X 5 


+ X 4 + X 2 





Z 


(To 

^ z 


00100 


00100 


1 




+x + 1 





(Tq 


(Tq 


01000 


01000 











(7/1 


(T/i 


10000 


10000 








^0 


(Tn -I- (T<? 


(Tn(Ti 


00011 


01001 








^0 


<7\ + (74 


(Ti (To 

^ 1 ^ z 


00110 


10010 








^0 


(To -|- (Tn 
^ z 1 ^ U 


(To (Tq 
^ z ^ 


01100 


00101 





X 5 


+ X 3 + 1 


^0 


(Tq 4- (Ti 


(Tq (T/i 


11000 


01010 











(T4 -|- (T2 


(T/i (Tn 

^ 4^ U 


10001 


10100 










(Tq (7 a 
w l ^4 


(Tn (To 


00101 


11000 











(T/1 ~h (Tn 


(Ti (Tq 


01010 


10001 






+ X 3 + X 2 





(Tn + (Ti 


(To (T/i 


10100 


00011 







+x + 1 


9 

On 




(Ti + (To 
^ 1 1 ^ z 


(TQ(Tn 


01001 


00110 










(To -4- (Tq 

z 1 w 


(T/i (Ti 


10010 


01100 








oZ 


fTrt -\- (To -4- (Tq -4- (T/i 
^ U l ^ z l w l ^ 4 


(Tn(Ti (To 
U 1 z 


00111 


11101 








(Tn 14 
^0 


(Ti -|- (Tq -|- (T/i -4- (Tn 
^ 1 1 ^ 1 ^ 4 1 ^ U 


(Ti (To(Tq 
f 1 ^ z^ 


01110 


11011 











(To + (Ta H - (Tn + (Ti 


(To (Tq(Ta 

^ z O 


11100 


10111 





x° 


+ x 2 + 1 


U 


O3 + Oo + Oi + 2 


0304 O"0 


11001 


01111 










<T 4 + (Ti + <7 2 + O-3 


O-4<T O"l 


10011 


11110 










°"l + °2 + 0"4 


a (Tia 3 


01011 


10110 








a 22 
°o 


2 + O3 + O"0 




10110 


01101 




X 5 


4 1 3 
+ X + X 


rx 13 


(J3 + (T4 + (7\ 


a 2 (T 3 a 


01101 


11010 


1 




+x 2 + 1 




(J4 + (To + 0"2 


o 3 o 4 Oi 


11010 


10101 








< 


C"0 + Cl + C"3 




10101 


01011 








rr 15 


O"0 + 0"3 + 0"4 


00010203 


01111 


11001 








30 
°0 


G\ + O4 + O 


01020304 


11110 


10011 




X 5 


4 1 3 
+ X + X 


rr 29 


0"2 + O"0 + °"l 


O"2O"3O"4O"0 


11101 


00111 


1 




+x + 1 


^o 27 


(J3 + <7i + (T2 


0-30-40-00"! 


11011 


01110 








-o 23 


0" 4 + 0"2 + O3 


O4CT0O1O2 


10111 


11100 








°o 


O + Ci + 2 + O3 + O4 


000^20304 


11111 


11111 


1 


x + 1 











00000 


00000 





X 



Table 6: 
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